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I present a detailed derivation of wideband optical pulses interacting with a Raman transition 
in the kind of scheme currently used to generate the ultra broadband light fields needed to create 
ultrashort pulses. In contrast to the usual approach using separate field envelopes for the pump, 
Stokes, and anti-Stokes spectral lines, I use a single field envelope. This requires the inclusion of few- 
cycle corrections to the pulse propagation. The single-field model makes fewer approximations and is 
mathematically (and hence computationally) simpler, although it does require greater computational 
resources to implement. The single-field theory reduces to the traditional multi-field one using 
appropriate approximations. 
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I. INTRODUCTION 

An important aim of current wideband Raman experiments is to try to efhciently generate few-cycle pulses 

If driven strongly enough, the two-photon Raman transition modulates the incoming field by adding sidebands 
separated by the transition frequency. Wideband fields are generated as these sidebands generate sidebands of their 
own (and so on), thus generating a wide comb of frequency components separated by the transition frequency. If 
a scheme can be implemented that adjusts the phases of each component appropriately, then few- or single- cycle 
optical pulses can be obtained (see e.g. Q). Standard theoretical treatments of this process split the field into fields 
components centred on the teeth of this comb. The approach has the advantage that the components can be modeled 
reasonably well with slowly varying envelopes, but of course it has the disadvantage of needing to keep track of a 
large number of components. 

In experiments like those of Sali et.al. the Raman transition is driven near-resonantly by a pair of intense 

pump pulses about lOOfs long; compared to the transition frequency of about 130THz, the spectra of each pump 
pulse (and hence the generated sidebands) are relatively narrow. This means that a multi-component model is 
still not unreasonable, even if numerical considerations might demand that the arrays used to store these spectra 
overlap in frequency space. However, if we were to move to shorter pump pulses, or to a single (much shorter) 
pump pulse with enough bandwidth to efficiently excite the transition, we would reach the regime where the "teeth" 
from the spectral comb significantly overlap. At this point, not only would we be forced to move from an SVEA 
(Slowly Varying Envelope Approximation) solution of the field propagation to a more accurate Generalized Few-cycle 
Envelope Approximation (GFEA) 0,13, but the utility of multiple field components becomes questionable. In this 
regime it can be advantageous to treat the field as a single unit, rather than splitting it into pieces. Note that this 
approach still differs from solutions of Maxwell's equations such as FDTD (finite difference time domain) or PSSD 
(pseudospectral spatial domain) [lO|, because our single-field is based on a second-order wave equation, and uses a 
convenient choice of carrier function to define a field envelope. 

Following these considerations, we now derive a single-field model for Raman generation, and, apart from that 
notable detail, follow an analogous path to that of Hickman, Paisner, and Bischel (HPB) ^|. In the model, we find 
that the coupling constants retain an oscillatory behaviour at the transition frequency, and that it is this that impresses 
the sideband modulation on the propagating field. Since the field is not only wideband, but contains significant 
sideband components, we need to propagate this (no longer slowly varying) field envelope using the GFEA. The 
necessity of allowing for these can be demonstrated by converting the single-field model into a multi-field counterpart 
- without the envelope-gradient corrections, we will not get a correct multi-field model. 

A. Summary of the theory and the numerical implementation 

We model the wideband Raman generation process in the following way. We specify the field frequencies {uji) of 
interest, which are usually at integer spacings of the transition frequency (llJa) from the main pump laser frequency 
(wo)- Each of these field components is described usin a standard envelope theory (i.e. as Ai{t)) with a time-history, 
allowing us to simulate pulses as they travel through the Raman medium. The Raman medium is modelled as a 



2 



two- level atom using a (its) Bloch vector (w, v,w), and this Bloch vector is driven by each combination of spectrally 
adjacent field components (~ ^ AiA^i). Each of the field components Ai is driven by the atomic polarization v) 
in combination with its pair of adjacent field components (i.e. Ai-i, Ai^i). 
Each field component evolves as 



d^Aj = — {-[«'- Aj+i exp (+^(fc' ^ - fc' )z - lAt) + [v' + tu'] exp (+i(fc'_i - fc' )z + lAt) | 



-{k,-ko)A, (1.1.1) 

The transition evolves as {p[2 = (u',v')) 

Vi2 = |-l2 + ^A + 2^g'Y,A*A^ p[, + 8^f'Y,A,A*_,.w.e+<''^-''^-^^^-^^' (1.1.2) 



dtw 



-71 (w - w,) + 4i/'^ [A*ylj+ip'i2-e+'^''^+'"''^'^^"''^* - A*Aj_ip'i*2.e+'('=^-i-'=^)^+*^* (1.1.3) 



Here a is the number density of the atoms or molecules; 71, 72 are the population and polarization decay rates for 
the transition; /' is the field-transition coupling constant; g' is the stark shift coefficient; coj^kj are the frequencies 
and wavevectors for the field components, A gives a rotating frame for pi2 — > 

Additionally, a Cauchy dispersion is applied to the propagation of the field components. 

Because each of the field components has a time-history, this translates to a spectral width. In typical cases with 
roughly nanosecond pulse lengths, the bandwidth of each component will be tiny compared to the transition frequency, 
so there will be large uneventful gaps in the total spectrum. In contrast, if the pulse lengths drop to (say) roughly 
100 femtoseconds, the bandwidths of the field components will form a noticeable fraction of the total spectrum 



B. A comment on Cauchy dispersion 



As regards the mismatch terms in the Raman model, Geoff New has remarked (email, 20040121) that there's a key 
point about the 7n's that is not made properly in most of the McDonald et al publications; since it is usually said 
rather enigmatically that the 7's are "parameterized" by 71. The point is that if one assumes a Cauchy- type law for 
the refractive index variation, all the gamma's are linked by a recurrence relation, and so you only need to specify 
one of them, from which all the others will follow. The point is made properly (to GN's knowledge) only in ref 0]. 



II. SINGLE-FIELD WIDEBAND RAMAN 



Note: This single field derivation does mean some of the approximations as to the "slowness" of the field variation seem 
somewhat stringent. However, since a conversion to a multi-field model ts possible, it would seem the field variation constraints 
are less stringent than they would first appear. 



A. Coupled wavefunction equations 



I start by considering the wave function V' of a single molecule (e.g. H2) and the electric field E, and write the 
time-dependent wave function by expanding it in terms of the eigenf unctions in the field-free (i.e. E — 0) case. 
This means I can get the expansion coefficients by solving for an effective Schroedinger equation that contains a 
two-photon Rabi frequency created by means of an interaction term based on a field-dependent dipole moment. I 
assume a dispersionless medium and write all equations in terms of position z and retarded times t = tiat — z j c. Here 
I follow the method of HPB ^Jli but use only a single E field rather than multiple components. Note that HPB use 
Gaussian units, so there can appear to be inconsistencies when comparing my formulae (in S.I.) to theirs. 
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I denote the known molecular eigenfunctions of the unperturbed Hamiltonian Hq as |n), and their corresponding 
energies HWn- I want to obtain the solution to 

{Ho + V)ij = ih^ (2.1.1) 

for V = -dE (2.1.2) 

V = ^c„e-^^"*|n), (2.1.3) 

n 

where d is the electronic dipole moment operator and the c„ are a set of complex probability amplitudes. 
A standard derivation for the equations of motion of the Cj co-efficients proceeds as - 



= {Ho + V)J2cje-^'^^'\j) (2.1.4) 

i 3 

^ft^|-^mc,e-^^* + e-^^*^||^) = (ifo + V) ^ c^e-'^^* |i) (2.1.5) 

th{n\ Y^LtWide-^"^^' + e-'"^''^] \i} = (n| {Ho + F) ^c.e"'^^* \j) (2.1.6) 

i j 

ihLzWnCne-''^-' + e-^'^'^*^\ = hCnWue-^"^-' + Y,Cje-^'^^' {n\- d.E\j) (2.1.7) 

^;i|-^W„c„ + ^| = hWnCn-Y.Cie-^^^'+^^-' {n\d.E\j) (2.1.8) 

3 

We now use perturbation theory, & dnm = {n-\d\n), following two independent (but related) strands. 

1. CASE (i): Electric field 

This strand leaves the perturbing potential as a function of electric field E, and does not replace it with a carrier- 
envelope description. Although apparently the simplest strategy, it is generally impractical as it imposes constraints 
on the field and other model parameters that are too restrictive to be useful. 

dc 

ih-^ = - Y,Cidr,iEeK:s)[-i{Wi-Wr,)t] (2.1.10) 

i 

... assume Cn^dni vary only slowly, so I can integrate just the exponentials (2.1.11) 

- -^}^c,rf.,i? (2.1.12) 
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Now substitute the c; solution into the dcn/dt equations, and introduce the shorthand notation Wij = Wi — Wj 



dcn 
l—r- 
dt 



E 



exp [—iWjit] 



W, - 



dniEeyip [-iWint] 
I x - x - eyiT)\—iWjit] , ^ r Tir 1 



p2 1 V- , , ^ exp [-iWjnt 



where 



i exp [-iWjnt] ^ 



(2.1.14) 
(2.1.15) 

(2.1.16) 
(2.1.17) 

(2.1.18) 



UNITS 



UNITS: [anj] 
a 



= CWJ-i =(jF-i)'mV-i - JmV-f2.1.19) 
Jm2y-2 J-ig^^ (Fm^i)^ = J.J^i X y^^.y^ X rn^.m^^ X s-i =5^^(2.1.20) 



,r\s-\Cm.Cm. [s 



Since we are only interested in the end states j = 1,2, between which the Raman transition occurs, we only need 
calculate ci, C2; however we still retain the sum over all i intermediates states, as they affect the coupling between 1 
and 2. The diagonal couplings {anj,n = j} are real; but the off-diagonal couplings {anj,n ^ j} undergo complex 
oscillations according to the difference in their energy levels. 

Their frequency dependence is discussed after the following subsection; HPB's corresponding parameters, which do 
not oscillate, were assumed to be frequency independent. 



2. CASE (ii): Electric Field Envelope 



This strand replaces the electric field E with a carrier-envelope description, but, unlike HPB, I use only a single 
carrier-envelope component rather than a set indexed by some integer j. This is necessary, because in the previous 
strand I ended up with coupling constants a„j with strong frequency dependences. 

I introduce the envelope and carrier for the field: 



so that 



ih- 



E 

dcn 
dt 



- exp [-1 {Wi - Wn) t] [Ae-'-""* + A*e+''^°'] 



(2.1.21) 
(2.1.22) 



now use / — > i; and assume c„, d„i vary only slowly, so I can integrate just the exponentials 



(2.1.23) 



, \ , exp [-1 {Wj -W,+uo)t] , ^ , exp [-i {W, - W, - u;^X t] 



24) 



1^ ^ r exp [-1 {Wj -W^+ ujq) t] exp [-1 {W, -W^- luq) C 



W^ - Wj + WO 



1.25) 



Note the swap of {Wj — Wi) to — (Wi — Wj) in the denominators. Now substitute the Ci solution into the dcn/dt 
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equations (using I j), and introduce the shorthand notation Wij — Wi — Wj, 
dcn 



ih 



dt 



where 



^^ exp [-^{Wj^ +uja)t] ^ exp [-i (Wj, - luq) t] 



Wij - LOo 



Wij - UJQ 



i 3 

, . [ .2 exp [-^ (Wjn + 2ujo) t] exp [-iWjnt] 

I .IM ^^P ^ exp [-Z (Ty,-„ - 2u;o) t] 



^^ exp +t^o)^] _^ ^* exp [-z (W'jn - luq) t] 



Wij + UJQ 



dmEex^p[-iWint] (2.1.26) 
d„, [Ae~'^"' + A* 6+'^"'] exp [-{flff^^S]) 

(2.1.28) 



Wij + ujo 



W^j + Wo 



1 



i ^ ^ CjdnidijAA* exp [-2Wj„t] 



W^j +uja + Wij - ujQ 
{Wij - too) [Wij + ujo) 



^ ^ CjdnidijAA* exp [-iWjnt] — 

i 3 

-AA* ^ Cj ^ jdnidij exp [-iWjnt 

3 i 



2W,j- 



i^nj = ^ exp [-iWjnt] 



= i exp E d„id; 



(2.1.29) 



by discarding the 2a;o terms; (2.1.30) 

(2.1.31) 

(2.1.32) 

(2.1.33) 
(2.1.34) 



(2.1.35) 



These redefined a'^j parameters still oscillate, as they must because unlike in the HPB derivation, there is no 
frequency difference between field components to cancel with the Raman transition frequency. The coupling also 
varies with frequency, which is discussed next. 



B. Raman coupling parameters: approximations 

I now discuss two particular (and vital) approximations applied to the Raman coupling parameters. 

First, note that (as in HPB), I will take the indices 1 and 2 to correspond to the two states involved in the (Raman) 
transition of interest; these will be the and 1 vibrational (or perhaps rotational) levels of the electronic ground state. 
Indices 3 and above will correspond to (quoting HPB) "translational motion on higher electronic states" . 

Note: I can see that assigning these indices to higher electronic states will conveniently keep the energy separations 
for transitions to greater than that of the 1 <-!■ 2 transitions, but it 's not so clear to me why I can ignore all the higher 
vibrational (or rotational) states. 

Since I am interested only in the Raman transition, I specialise the above equations for the coefficients c„ , calculating 
ci and C2 only, and assuming that the di2 = (11^1) dipole moment is zero. This means we will only be including 
transitions between indices 1 and 2 that go via one of the higher states j > 3, since we still allow dij, d2j ^ Q Vj > 3. 
Further, I solve for the coefficients for the higher states in terms of ci and C2, in an adiabatic approximation justified 
when ci and C2 vary only slowly compared to the exponential terms. 

Note: Separate from the oscillations that occur in my coupling parameters (but not in HPB), there is the issue of 
frequency dependence which applies to both HPB and my parameters. 

For both HPB (their eqn.(15)), and my field-carrier based derivation fean. (|2.1.35|l ). the presence of the field carrier 
in the denominator is helpful. Since it can reasonably be assumed to be much greater than the inter (Raman)-level 
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energy differences, the fractional terms will be correspondingly small, and so the a„j (a'y ) parameters can be assumed 
to be independent of frequency (or at least nearly so). 

This is not in any way obviously (or even partially) true for the no-carrier E based fean. (|2.1.18|l . Unj) parameters, 
which depend only on the difference in Raman levels - at best you might perhaps hope that the denominators were 
approximate multiples of each other. 



1. Near-zero difference 



Here I assume that the forward and backward transitions (between levels 2 and 1) have nearly the same amplitude. 
It is not self-evidently true, but HPB must have made an equivalent assumption. In any case, since — is the 
difference of similar terms, it will be smaller by at worst a factor of two ((1 + (5 — !)/(! + + 1) ~ S/2) 



'-12 



4i = Texp[+iW2it]'^dlid*. 



W,, 



i2 



i2 



i2 



- i exp [-zM^i2t] ^ d2idix 



il 



exp [+zW2it] 



(2.2.1) 

(2.2.2) 
(2.2.3) 



This approximation allows me to equate a'jj^ to a'12, and hence change — a2]^c*ci + a'*2C2C2 — > a'21 {c^C2 
the p\2 equation below. This is a vital step in getting to a simple form equivalent to the Bloch equations. 



c{c\) in 



2. Sum is double 



This follows from the above assumption (as per HPB) that the forward and backward transitions (between levels 
2 and 1) have nearly the same amplitude. This approximation allows me to replace a'21 with a[2 in the w equation 
below, which simplifies my notation and makes the coupling term match that in the pi2 equation, important in getting 
to true Bloch equations. 



"12 + "2*1 = exp [-iW2it] ^ dud: 

i 

, , r w,2 



W,2 
»2 Ti/2 , ,2 



i exp [+iWi2t] d*2,d*i 



Wf2-^l 



exp [-iW2it] 



Ot\2& 



(2.2.4) 

(2.2.5) 
(2.2.6) 



3. Sum is cosine 



I do not reduce the exponential sum to a (trig) cosine function, for potential use later when matching and/or 
differencing exponentials. However, this sum is does not occur, so the 012 is never used. 



1 = i exp [-iW2it] ^ dud: 



»2 Ti^2 , ,2 
^i2 - ^0 



exp \-%Wx2t\ ^ d2idii_ 



T dudi2 



h 

0.12 |e 



■ exp [-zVF2it] 



' 12-^0 



Wi2-^0 



■ exp [-l-iVF2it] 



(2.2.7) 

(2.2.8) 
(2.2.9) 
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4- Stark Shifts 



Here I calculate the Stark shift term 



= ^ cxp [-iWiit] diidg 



-y 



W,2 



i exp [+lW22t] d*2^d, 



w, 



i2 



i2 TT72 , ,2 



Simplify assuming W21 ^ W^ii, M^i2, while assuming W^i/ {Wf^ — Wq) 1; or, more accurately, that 

W^2"l 



(2.2.10) 
(2.2.11) 



1 > 



for n 



1, 2, 



and keeping only the terms first order in W2i/Wii, W2i/Wi2, etc 

WiA . , ,2 Wil + W21 



-y 



M2. 

|d2. 
\d2i\ 



WA+2WaW2i + Wi^-uj^ 
2 W^i + W21 



{Wl - ul) [1 + {2W2,/W,, + WiJWl) Wll {Wl - 



,2 + W21 

2 Wg + W21 



1- 2 



1 - 2 



.^21 W^li 
W21 wl 



\dli\^ - \d2^\^ - \d2i\^ 



2 W21 



2\d 



2i\ 



(2.2.12) 

(2.2.13) 
>2.2.14) 
(2.2.15) 
(2.2.16) 

(2.2.17) 



Ignoring the transition frequency terms, this is just the difference in the energy shifts of the two levels with 
field intensity. This is a purely real quantity, with no imaginary part. 



C. Two photon Bloch equations 



From now on, I use 3 (a script "I") to represent if following on from III Ka\ and eons. H2. 1.171 I2.1.18|l : or to 
represent 2 \ A\ if following on from III A2l and eons. (|2. 1.341 12.1.!?^ . Since it is the envelope-carrier description of the 
field E which is most useful fin III A2ll . for most purposes [J = 2 |A| holds; remember that it is difficult to maintain 
the accuracy of the approximations relied on above for the E"^ picture fin III A l]) . I also drop the prime on the a'^^ 

parameters used in the (Case lii)) "2 |^|^" electric field envelope model. 

Since we are only interested in ci , C2 , and because we keep only stationary or slowly varying terms, we can write 
equations for ci , C2 as - 



dci 

It 
d4 

~dt 

dc2 
dt 

dc* 
in — - 
dt 



-aii3ci — q;i2Jc2 
-a2i3ci — q;22Jc2 



In matrix form: 



d 


Cl 


1 


and a.12'} 




Cl 


di 


C2 


ih 


a2i3 a223 




C2 



(2.3.1) 
(2.3.2) 
(2.3.3) 
(2.3.4) 

(2.3.5) 
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D. Alternative formulations: Kien et.al. (KLKOHS), Hickman et.al. (HPB) 



This interlude compares my matrix equation for ci, C2 to those from Kien et.al. [13 (KLKOHS) and Hickman et.al. 
[ni (HPB). This is useful as a fixed point of refernce beween the approaches, enabling quick conversions between the 
parameter variables. 



cf. (KLKOHS) eqn.(13): 
where Bq are KLKOHS's Eq envelopes: E 

n.n. 



d 


Ca 


1 


h^^aa fl^ab 




Ca 


di 


Cb 


ih 


Klba fi^bb - S 




Cb 



Y.^,BqB*^ 



h^ab - -J2dgBqB*q = 

fi-^bb = - bqBqB*^^ 



ha 



(2.4.1) 

(2.4.2) 
(2.4.3) 

(2.4.4) 

(2.4.5) 



cf. (HPB) eqn.(13): 
where auij are HPB's parameters: 



d 


Cl 


1 


'-Hn 


—H12 




Cl 


di 


C2 




—H21 


—H22 




. . 



H12 = 
H22 = 



j 
j 

-\Y.^H22{iO,)V,V; 



- ^21 



(2.4.6) 
(2.4.7) 
(2.4.8) 
(2.4.9) 



Thus comparing the Rabi-like parts of my eqn (|2.3.5|l with that of KLKOHS (my eqn (|2.4.1|l ') and that of HPB (my 
eqn H2.4.6|l ). gives 



ai2'i = ai2.2A*A^hf.2A*A = | ^ - d,. 



--Y,aHl2{L0J)VJ^^;:}A■n 



• This uses KLKOHS: E=^\{Bi + Bl +B2+ B^) 

E^ ^ \{Bf + Bf + Bl + Bf + 2BiBl + 252^^ + 2BiB* + 25*^2) 

£/l_2 ~ 4^i3li32 — 2-°l-°2 

which tells us the size of the field contribution to the KLKOHS Rabi-like term. 

• This uses HPB: which tells us the size of the field contribution to the HPB Rabi-like term, (cf KLKOHS). 

So U,aa — o:ii3/h, Qab = ai2J/ft, ^bb = a223/^ (temporarily ignoring their detuning d); and looking ahead to the 
definition of / = ai2/2h = {ai2 + a2i) /2h gives us f^afc — /3, noting d — 2AqA*_^j^, compared to 3 = ^EqE*^-^^. Note 
that KLKOHS have defined ^ab = 5 dqEqE*^^^, so that f = dq. 



E. Two photon Bloch equations (cont) 



An important difference between KLKOHS & HPB, and my equations is that my ofy coupling parameters retain a 
'slow" time dependence at the Raman frequency W12. I now turn my equations H2.3.5|) into Bloch equations by first 
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working out d/ dt of cjci, C2C2, and c\c2, 

ih — = ~T~ + — -i — — Q!ii!}c*Ci — Q!i2Jc*C2 + a*]^Jcic^ + a*2JciC2 (2.5.1 

(JjL LLv (JiL 

= +a*2JciC2 — Q;i2lIc*C2 — [ail — ctii] 'ic\ci (2-5.2 

= +3 (q;i2CiC2 — Q!i2CiC2) ; since an = aii (2.5.3 

^^ ^'^dt^ = lhC2^^ + ^^'^2-^ = -a2lJc2Ci - a22Jc2C2 + a2iJc2cJ + a223c2C2 (2.5.4 

— — a2iJc2Ci + a2i3c2Ci — [a22 — 022] 3C2C2 (2.5.5 
= — 3 (a2iC2Ci — a2iC2Ci) ; since a22 = 022 (2.5.6 

t?l — -^-^ = ihci—r^ + lhc2—T^ = +a2iJciCi + a22JciC2 — aiiJCjCi — ai2Uc2C2 (2.5.7 

Cifc LLL LLL 

— — (an — a22) 3c2Ci + aji^ciCi — ai23c2C2 (2.5.8 
= — (an — 022) 3ciC2 — ai23 [C2C2 — cjci] ; since a2i ~ ai2 (2.5.9 

Use P12 = C1C2 and w = C2C2 — CiCi, 

ift— = — J (a2iC2Ci — a2iC2Ci) — J (ai2CiC2 — ai2CiC2) (2.5.10 

= 3 (— aj2CiC2 + a2icjc2 — a2iCiC2 + ai2CiC2) (2.5.11 

= — 3 (aj2 + a2i) C1C2 + 3 (a2i + ai2) c^C2 (2.5.12 

= -3(ai2 + a2i)/3i2 +3(a2i +ai2)pt2 (2.5.13 

ih = — (an — a22) 3ciC2 — ai23 [C2C2 — c*ci] (2.5.14 

= - (an - a22) 3pi2 - ai23w (2.5.15 



Hence 



dpi2 (ail-a22)T , "12 a /orciR 

— — = I Jpi2 + z— — (2.5.16 

at n n 

dw _^ (a|2+a2l) (021+012) r,* /or; 17 

— = +z 3pl2-^ 3pi2, (2.5.17 

= +^Ml3p,2-^'^-^'ipl2■ (2.5.18 



F. Transformations of the Bloch Equations 

I define a new coupling parameter /' (c.f. ws), following the definition of ai2 in ean. (|2.2.6|l . and cjb = W2 — Wi = 
W2i- I also define a wa, which corresponds to the intensity dependent shift detuning shift proportional to (a22 — an). 
Thus, factoring any complex phase of ai2 into the angle 5' , I get - 

~ NB: this is a real quantity, see eons. 12.2.111 Ti. 2. 1711 (2.6.2) 

n 



9' 



^ («12+ oSl) j ^ ^^-^.,t-^5'J ^ ^-^u,t-^8' j (2.6.3) 

= (^il^^3. (2.6.4) 



The Bloch equations (|2.5.16l I2.5.17|l can now be rewritten 



^ = zg>i2 + i/W— (2.6.5) 
dt 



^12 6+*'^"*+'*' - ^^2 6"'"""*"'*' 



(2.6.6) 
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In analogy to both the standard two level atom Bloch equations, and those in HPB (HPB 18), these equations can 
have losses introduced, in - 



'^^'^ = -12P12+ 19^12 + if "iwe-^"^^'-''' (2.6.7) 
= -71 (w - w,) + 2if'3 \p,2e+"^'*+''' - ptse-"""*-*^'] . (2.6.8) 



dt 
di 
dt 



In the standard atom-field case, the atom and field frequencies are similar, so an atom carrier could be chosen to 
match its evolution the field carrier frequency, leading naturally to a detuning term. However, in the Raman situation, 
the atomic frequency is far removed from the field frequency, so the frequency evolution that appears in the definition 
of?/' f eon. H2. 1.31) ) is sufficient. 

Here I define two things (a) new u,v variables that represent the density matrix element pi2, and (b) allow for a 
"detuning" rotation in pi2. This detuning rotation looks rather like a carrier+envelope representation for pi2, but 
without a "-I-C.C." since pi2 is a complex quantity. There is no need for any kind of carrier+envelope transformation 
for the inversion w. 

So, 

I now adapt the Bloch equations to allow for the rotation in pi2, and introduce a detuned transition frequency 
uj'^^ujb- A, 

Ot ( Pi2e I = -«Api2e + e dtPi2 (2.6.10) 

= -72Pi2e"''^*"'*' + ig'Vi2e"''^*""'' + if"iwe-"^^'~'^' (2.6.11) 
dtw = --ii{w-Wi) + 2if'3(^p'^2^-'^'^'^'(i'^"'^''^"^' - 912^'^'^'^"^'^^^^^ (2.6.12) 
=^ e-'^'-'^' dtp\2 = {--f2 + iA)p\2e-'^'-'^' +ig':]p[2e^"^'''^' +if"Jwe-"^'''-'^' (2.6.13) 

dtW = —fiiw-w,) + 2tfd(^p[2e+<'-p[*2e-'^''') (2.6.14) 

=> dtp\2 = (-72 + *A) p'i2 + *.9'Vi2 + ^/'^e'"^* (2-6.15) 

dtW - -71 (u' - mO + 2i/'J (p'i2e+'"^* - P'i*2e~^"^*) (2-6.16) 

Notice that we have made the fixed complex phase (the 6') vanish from the equations; this is not dependent on the 
presence of a finite "detuning" A. In what follows, A = 0, and there are two cases depending on the chosen meaning 
for J; however note that both give the same result. 

CASE (i): 3 = E^: continue by applying the carrier+envelope for the E field E = Ae"'""* + ^*e+'"o*. 

CASE (ii): 3 2A*A: we can jump straight to ean. (|2.6.25l I2.6.26|l (i.e. (£>)), since the RWA was made when 
deriving the couplings anj- 
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The CASE (i) 3 ^ E'^ starting point is - 

dtp'i2 = i-l2 + + ig'3) +if'3we-<* (2.6.17) 

dtw = -71 {w - w,) + 2i,n (p'i2e+""^* - p'lV^'') , (2.6.18) 

iB)=^ dtp[2 = {-l2 + tA + ig'E^) ifE'^e^'^Kw (2.6.19) 

dtw = -71 {w - w,) + 2iJ'E^ (p'i2e+^"^* - p'i*2e"*"^*) , (2.6.20) 

{C)=^ dtp'^2 = (-72 + «A + p'i2 +if'[Ae-"^°' + A*e+''^°']'we-'<' (2.6.21) 
dtw = -71 (w - !«,) + 2if' [Ae-^^o* + ^*e+''^«*] ^ (pi2e+*'^^* - pt2e"*'^^(^)Q-22) 

(apply a RWA about Wo) [^e"'""* + A*e+''^°*] ^ 

= A2e-2*-'«* + 2^M + ^*^e+2*"«*l (2.6.23) 

=^ « 2^M (2.6.24) 

C^S'£;(m) & (D) =^ 9tp'i2 = (-72 + «A + 22.9'A*A)p'i2 + 2i/'A*Awe-"^^* (2.6.25) 

- -71 (?«-«;,) + W'AM (pi2e+^'"^* - pJae"'"^*) , (2.6.26) 

(£;) 5tp'i2 = (-72 + iA + 225'^M)p'i2 + 2«/'v4*Aw.e-'"'* (2.6.27) 

= -71 (z« - i«.) + ^if'A*A (p'i2e+'"^* - P'i*2e"''"^*) , (2-6.28) 

(F: split pi2,)=> dtu' = -72W' - (A + 25'A* A) i>' + A* Aw.fiin{uj'^t) (2.6.29) 

9ti;' = -72i;' + (A + 25'A*A)u' + 4/'AMt(;. cos (w[t) (2.6.30) 

dtW = --ii{w'-Wi) -A,f'A*Au'.iiin{uj't^t) - Af A* Av' .coii{uj'^t) , (2.6.31) 

Note 2pi2 = u + IV, since pi2 = cicj, as per ean. H2. 6.9(1 : this means the 2 in the m, u equations becomes 4, whereas 
the 4 in the w equation isn't doubled to 8. Note (again) that the Stark shift parameter g' is real valued. 

1. Rotations 

Looking at eans. H2.6.29l I2.6.'5ni 12.6.311) . we can see three rotation angles: Ouv — A, Ouw — Af A* Asm{uji,t — 5'), 
and Ovw = Af A* Acos{uJbt — 5'), which apply to the coordinate pairs {u',v'), {u',w), and {v',w) respectively. For 
extra generality, I will allow a complex valued f = fr + ^fi, but note that previous definitions make /' real. ) 

In vector notation, ignoring the losses, we can construct a torque vector O (unchecked signs), 

dtu' = O.u' - Ouv-v' + euww (2.6.32) 

dtv' = A.u' + 0.v' -9yn,-w (2.6.33) 

dtw = -Ouw-u' + 0yn,-v' + 0.W (2.6.34) 

^A^'^v',w] = [-AflA*A, -AflA*A, A] (2.6.35) 



dt 



= [Ovw, Ouw, 6'„„] X [u', w', w] (2.6.36) 
=^ = (ixU (2.6.37) 

I should now be able to turn this into a (unitary) rotation matrix for the U vector. 

G. The polarization driving the field 

There is a distinction between the atomic polarization P, and the effect on the field of that atomic polarization. 
This is beause we are dealing with a nonlinear interaction. I denote the (effective) Raman polarization T, and this 
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quantity is the one that appears in the wave equation. AUen and Eberly "Optical resonance and two level atoms" [l5l | 
have, for the standard (non- Raman) case (skipping the integral), an eqn.(AE 4.2) - 

(^£■4.2) P{t,z) = m[ucos{Lut- Kz) -vsm{ujt- Kz)] (2.7.1) 

HPB, after summing the electric field components (which absorbs a 1/2), have 

{HPB 21) P = ie-*^(ii + iv)(7ai2^ [V,e*'^^-i* + V,*e-*'^^+i*] + c.c. (2.7.2) 

3 

= ie-'^(ii + iv)aai2e-''^'*i^ [V,-e*'^'* + y/e-*'^^*] + c.c. (2.7.3) 
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= ^e-''^ {u + iv) aauC"^"* E + c.c. (2.7.4) 

I now write down the the polarization envelope i?, in my variables, and based on the same carrier as the electric 
field. I also introduce a complex factor equivalent to HPB's e~*^, for closer matching of the expressions. 

y ='Be-'~ + 'B*e+'~ = piae^'Vaia (Ae^*- + A*e+'") + c.c. (2.7.5) 
= Cpi2e"*V.ai2e~*'^^* (yle-*- + A*e+'-) + c.c. (2.7.6) 
= Cpi2e"*V.a'i2e~*'^^* {Ae-'~ + A*e+'-) + C,p\*c^-'^a.c!^^c^''<^ (A*e+'- + ^ffTf^) 
= Cpi2e"'V.a'i2e+''^^*^e-*^ + Cp'i2*e+'V.a'i2e"*"^*Ae-*^ 

+Cp'i2e~'^cr.a;2e+'"^*^*e+'- + Cp'i2*e+'V.a'i2e""^^*A*e+'= (2.7.8) 
= Cc^-ai2|Pi2e " +Pi2 e " ^ | 

+CCT.a'i2 {^'126+'"^*"'" + p'i2*e"*'^^*+"'} A*e+^^ (2.7.9) 
= CCT.a'ia^W^e-'" + Ccr.a'i2^W^*e+*^ (2.7.10) 
where the real valued is X{t) = i^p[^e+''^'''*-''^ + ^'12*6^*"^*^'"} (2.7.11) 

= l^^{u' + iv')e+"^'^*-''' + ^{u' ~iv')e-"^''''+'''^ (2.7.12) 

= {u' cos {io't,t^0)-v' sin {uj't,t~0)} (2.7.13) 
and 'B{t) = (aa[^X{t)A{t) (2.7.14) 

From the post-o^o RWA at eons. (|2. 6. 291 12.6.^ l2.().3Tll : above; and an un-scaled eqn. 

dtu' ^ -j2u' ' (A + 2g' A* A) v' +4:f'A*Awsm{uj'i,t) (2.7.15) 

dtv' = -"f2V + {A + 2g' A* A) u' + Af A* Aw cosiuj^t) (2.7.16) 

dtw = -7i(w'-w,) - 4/'AM. sin(w^t)u' - 4/'^* A. cos(w^i)w', (2.7.17) 

d,A{t) = NB: units [A]/m = ^^[03] ^ [A] = ^ [x] = [S-^] = 1 (2.7.18) 

Co no 47reo m.s ^ 

23(0 (2.7.19) 



2co7ioeo 

Co-ai2t^o 

I— 

2conoeo 

Co-ai2a;o 



[u cos (cjfct - 61) - w' sin (cjfct - e)] . A{t) (2.7.20) 

- 2conoeo ' ^^^^^^^^ 
Comparing the prefactor of ean. H2.7.21|l to that of the corrected (HPB 22), 

(HPB 22) 5^ = fe-*« {u - iv) V.+i ~ e+*« (u + iv) (2.7.22) 

dz 4eoc L ^ / J J 

we can see that the only apparent differences are a factor of no, and in that HPB have a e^*^ term that I omit. By 
imagining the cosine term split up into + and — frequency exponentials, it is easy to see how the relations between 
Vj and Vj±i will arise. Note the appearance (in HPB) of a carrier-dependent uji term, which in a standard multi-field 
approach would lead to different prefactors on the different d^Ai equations. 
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1. Simulation "photon" variables 
Scale the field E into square-root intensity "photon" variables, from / = 2noCoeoE'^ , so that 



Hence fA*A 



since 



Ap{t) = V2conoeo • Mt) 

r 



2conoeo ' 2corioeo ' 

fpApAp, so that 

-72u' - (A + 2gpA*A) v' + 4:fpA*ApW sm{uj'^t) 
-72W + (A + 2gpA*A) u' + AfpApApW cos{uj',^t) 
-71 (w' - w,) - 4:fpA*pApu' sm{uj'i,t) - AfpA*pApv' cos{uj'i,t) 
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dtu' 
dtv' 
dtw 



d.Ap{t) 



-Ap{t)X{t) 



2conoeo 

4fecnnoeo/p 



2conoeo 



Ap{t)X{t) 



2^{Aah)LUofpAp{t)X{t) 



ai2 



2hf 



ai2 



2corioeo 



2hfp 



i^uJoRfpAp{t)X{t); 



ai2 = AhcQUQeofp 



R = Aah 



5.31 X 10" 



Scahngs: 

h = 1.05 X 10-34 Js -> 1.05 X 10-i°nJ.fs 

a = ^ X X lO-iS/im-3 

2conoeo = 5.31 x lO^^^^ g-i j.m-i.V^^ 

Scalings: 

[ceo] = TO-s"^ ■ = J.S-1.C2.J-2 

10-24.C2.nJ-i./s-i 



(2.7.23) 
(2.7.24) 

(2.7.25) 
(2.7.26) 
(2.7.27) 

(2.7.28) 
(2.7.29) 
(2.7.30) 

(2.7.31) 



J.s-i.V-2 ^ 5.31 X 10-9 . nJ.fs-i.V-2 

= C2 . lO-^nJ-i . 10-15/s-i 
lO-^nJ-i . 10-V"^"^ = lO-i^C^.nJ-i 



./im 



Comments on the field propagation 

We can see that the derivative of A is proportional to lA, hence the (time-domain) evolution just rotates each point 
A(t) differently without changing its amplitude. This might seem to imply that we will never get a shorter pulse than 
we put in; but note that the A is an envelope, and the field is A -I- A* , so that in principle the phase of A might be such 
that its amplitudes cancel in certain t regions but not others, leading to a shorter pulse (and the magnetic field H 
also, since the fields are plane polarized). However, it seems unlikely that this will happen (barring some miraculous 
coincidence) from purely Raman effects. 

If I can predict the output spectral phases though, a structure with a suitably designed dispersion might be able to 
impose the desired phases. It would only be necessary to get the dispersion right at the comb points of the spectrum. 
See Shverdin et.al. jlfi] . who do a four-wave mixing optimization procedure in their experiment to match their phases 
appropriately. 



H. The steady state and gain co-efficient 



Starting from the eqns. assume v' is steady state (NB 72 — I/T2), so (with w 

dt~Q = -72w' + 4/'A*A(-l)cos(tJ6i) 



4r 

72 



A* A. cos{ijJi,t) 



4/'T2 . A* A. cos{LUbt) 



-1)- 

(2.8.1) 

(2.8.2) 
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I now use the field propagation equation H2.7.30() . with the A field envelope on each side scaled into "photon" 
variables. Inserting the steady state of v' calculated immediately above, we have - 

dzAp — i^^^^^-^^Apv'n sm{u}bt) (2.8.3) 
2conoeo 

= /'^°'i2^o A 4j'2^2 . A*A.cosiujbt).sm{LJbt) (2.8.4) 
2conoeo 

= ^^-P^Ap . m^:^ [cos{u;,t).sin{c.,t)] (2.8.5) 

- , ApA;Ap. [cosiiobt). siniiobt)] (2.8.6) 



= iGpApA*Ap.[cos{uJtt).sm{ujtt)], (2.8.7) 

Apparently, therefore, this Gp corresponds to the usual "gain co-efficient" g: 

cf(SMN) g'^^-^f^I^, (2.8.9) 

and we can assume that they are identical but for a factor of In^/C,. I do not know whether SMN silently assumes 
no = 1; or factors no into c. In any case no = 1 is usually accurate enough in gases. 

However, one hidden complication with my above equation is that Ap includes both the center "pump" field and 
the Raman sideband we are amplifying - thus to properly check the gain co-efficient, we should expand it into its 
components, as is done below. 

1. Pump and sideband calculation 

Ap includes both fundamental Ai and its Raman sideband A2 — For the moment, I leave the calculation 

in a rather abbreviated state - probably it should be shifted to the Multi-field section following later. So 

ApA*pAp = [Ai + A2] K + A^] [Ai + A2] (2.8.10) 
= [AiA{ + AiAl + A2AI + A2A*2][Ai + A2\ (2.8.11) 
= + AiA^Ai + A2AlAi + A2AIA1 + ArA\A2 + AiA^Aa + A2A\A2 + ^K'^^ 

1st order terms in A2 only: ~ ^1^2 + '2A2AIA1 (2.8.13) 

drop counter rotating AjAj: ~ 2A2AIA1. (2.8.14) 



Thus the sideband envelope A2 evolves as 



dzA'2 = iGpA'2.2AlAi.cos{LJbt) (2.8.15) 
= iGpA'2.Ii.cos{tJbt) (2.8.16) 
= iG' A'2.Ii.cos{uJbt), (2.8.17) 



since Im = h^m^m- I don't convert A'2 into I2 because it occurs equally on both side of the equation 



Gp (2.8.18) 

;Q/tQtQ/- 



III. MULTI-FIELD VARIANT OF SINGLE-FIELD RAMAN THEORY 

The single-field Raman model above can be converted into a traditional multi-field model as developed in e.g. HPB 
[lll | or Syed, McDonald and New |0| by replacing the field envelope with a sum of multiple envelopes using carrier 
exponentials spaced at the Raman frequency. When doing this, I will only get the correct multi-field form if few-cycle 
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(either SEWA or GFEA) corrections to the field evolution part of the theory are applied to the effective polarization 
caused by the Raman transition. 

The idea is to replace the single field envelope with a sum of multiple envelopes spaced at the Raman frequency, 
which are best placed to represent the comb of frequencies generated by the Raman interaction. Note that it will not 
necessarily be identical to HPB and/or SMN, because the field equations are derived from a propagation equation 
using a LUo,f3a carrier, but it will be very closely related. 

Starting from eons. (|2. 6.2712.6.2^11 . Since the single-field evolution equation fean. (|2.7.21|l ) uses an envelope A that 
is based on a carrier (see eon. (|2. 1.211) '). the single- field envelope A is replaced with Aj^s at frequency luj — uiq + juJt 
with wavevector kj = k(u)j); thus lo'^ = ujj — ujq, kj — kj — k{LOo) = kj — fco! also f3 k. The single-field envelope in 
terms of the new A, 's is 



A = Aj exp [—1 {u!jt — kjz)] 



(3.0.1) 



A. Polarization (p'12) 



First I will handle the polarization {P12) equation (from ean. H2.6.25() ^ (watch for any uj'f^ — uji, — A confusion, and 
note 2p[2 — u' + iv') - 



{Al) 



{RWA) 



dtp'^2 = (-72 + «A + 2ig'A*A) p'i2 + 2if A* Awe-^"^'^* 
dtp\2 = (-72 + *A)p'i2 + 2*ffyi2EE^l^'=^'''^"^""'^*e+^^-'=^+'='=)^ 



2 k 



-72 + lA) p'12 + 2w'p'i2 E E A*Afce+'(^-'=)"''*e+*(-'=^+'='=)^ 



2*/'EE^;^'»^ 



^+i{j — l — k)uJi,t+iAt+i{kk~kj)z 



(3.1.1) 



(3.1.2) 



(3.1.3) 



j k 



(-72 + *A) p',2 + 2^g'p'l2 E ^1^^"' + E A*A,^,.w.e+^^Ke+<''^^--'^^^^ (3.1.4) 

3 3 

72 +^A + 2^5'E^J*^J ) p'12 + W E ^J'^J+i-^-^^'^*-^"^'''''"'''"'^' (^.1.5) 
j J j 

72 + ^A + 215' E AjA^j P'i2 + Ww.Me E ^j-4*+i-e+'^*.e+*('=^-'=^-i)" 

E AjA*+^.e+'^*.e+''^''^-''^-'^' 



(split ^12) dtu = -72W- I^A + 2.g'E^i^j ) v-8fwlm 

dtv = -72W + A + 2.g' E A*jA. 



tAt Q-i-^ikj -~kj — i)z 



(3.1.6) 
(3.1.7) 
(3.1.8) 



Where the factor 2if' turns into 4i/' because the double summation gives two identical terms that only occur 
once in the single summation. When split into equations for u' and v' , the corresponding factor becomes 8if'. This 
equation for pi2 is equivalent to (SMN 2)|l7l|. except I have just kj — kj^i whereas they have Aj = fcj — fcj_i — fco + fc-i; 
however note they have the reverse sign in their carrier wave, so the only physical difference is the fco — part; also 
my definitions of the coupling differs slightly. 
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B. Inversion (w) 



And now the inversion (w) equation (also from ean. (|2.6.25ll ) (watch for any uJb vs lu'^^ confusion) 



dtw - -7i(u;-w,)+4i/'A*a[p'i 



Pi2e ' - pT2e-'"^* 



(3.2.1) 
(3.2.2) 



(A2) dtw = -7i(w-mO+4i/'^^A*Afc.e+*("^-"'=)*.e+*(-'=^+'='=)^ 

j k 

= -71 [w - wO + 4i/' ^ [A*Afcp;2e+'^^"''+^^'^''*"'^* - yl*Afep;*2e^*^^"''"^^'^''*+'^*] .e+*(''^"t-3'-^.3) 
J" fe 

(WA) « -7i(w;-w,) + 8z/'^ [A*Aj+ip;2-e~*'^*.e+^('=^+^-^'^)"- A*Aj_ip;Ve+''^*-e+'^''^-'"''^'^"] (3.2.4) 
= -71 (u; - mO + 4i/'^ [^^^j+i ("' + .g-'^^.e+'^'^^+i-'^^)^ 



-71 (u; - m,) + 4iw7'^ [yl*Aj+i.e~'^*.e+^('=^+i-'^'^)^ - c.c. 



4w 



+ C.C. 



-71 (ly - - Su'/'Elm [A*Aj+i.e-*^*.e+'('^^+i-'=^ 



)2 



8i; 



(3.2.5) 



(3.2.6) 



(3.2.7) 



Note the RWA's above (for both polarization and inversion equations) discard modulations at frequency of multiples 
of LOh- Quite a lot of physics has been removed by these RWA approximations, although it is very reasonable except in 
the very wideband limit. For example, the effect of next-nearest neighbour field components acting on the transition 
have been ignored, as have all more distant field-field interactions. In the next-nearest neighbour case, the dropped 
terms would impose a rapid LOh oscillation onto the polarization pi2, which would in turn tend to impose sidebands 
at onto each field component. It is reasonable to ignore such sidebands in the narrowband limit studied by 
most users of a multi-field Raman theory; but, in principle one might extend a multi- field theory to include them by 
inventing a scheme to apply the sidebands to the field component they are (near) resonant with. 



C. Fields (Aj) 

Note that the field evolution equation already has a carrier of exp \—i {ujot — kQz)] factored out of it. Thus I use A' 
from A — A' exp [— i {ujot — koz)], not A. Finally, I need to insert the GFEA few-cycle correction to the polarization 
term, because my (j ^ 0) sub-envelopes Aj have an ijuht time dependence that cannot be neglected. 

From eans. (|2. 7.211 12. 7.11|l . I get 



d.A'it) 
dz E 6xp [—1 (uj'^t — k'j 



z = I 



1 + ^ 

LUq 

1 + ^ 

LUq 



^-^^A'{t)X{t) 
2eoCo 

Cc^^o«i2 
2eoCo 



(3.3.1) 



p'i2e+^"'* + p'i*2e"''"^* 



E exp [-« {ujjt - k'^z'j^.3.2) 



(match ujj terms) [ikjAj + dzAj] e = 



'^"^"12 {A,+ip[^ exp [-^^(fc^+l - k'j)z - lAt] 



2eoCo 



+Aj^ip'*2 exp (+«(fcj_i - kj)z + lAt) } 
X [ojQ + idt] e-'^'^"* (3.3.3) 
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Then, using 

[uJo + idt]e-'''^»' — > t^oe"'^""* - zVbe"'^'^''* — > [t^o + J^h] e-'^'^^* w^-e"*^'^''* (3.3.4) 

So 

[ik)Aj + d,Aj\ e-'^'"'* = {^j+ipia exp [+t{k'^+i - - lAt] 

ik'jAj + d^Aj = ^'^^^^u [Aj+i exp - k'^)z - lAt] + Aj^i exp [+^(fc^■_l - /s^)^ + lAt]] 

- ^^p^y ^Aj+i exp (+i(fc^+i - fc^)z - lAt) - Aj^i exp (+i(fc^-_i - k'j)z + lAt)] (3 

dzAj = |_ j„ _ exp (+i(fc'ii — k'j)z - lAt) + [v + m] exp — k'Az + i 

4eoCo \ J I \ J J 

— I (kj — fco) (3 

This is in agreement with both HPBjlJJ barring the opposite sign on the RHS - similar agreement occurs with 
SMN[i3 once I identify q — [v + lu). Note that generally A = 0, as it just controls a frame rotation for pi2. 

Note that we can assume, quite reasonably, that the multi-field envelopes Aj will be better behaved than the 
single-field envelope A. However, we have made more approximations, notably by RWA'ing away all the off-resonant 
cross terms driving the atomic transition so a multi-field approach is not always better. These off-resonant terms are 
at 2loo ± - see just prior to the starting point above of eans. H2.6.25l l2.6.2T)jl . 

In photon variables, the above field propagation equation is (using R = Aha, also note ^ = 2 to conserve energy) 

dzApj = -Rcjjfp {-[v - lu] Apj+i exp {+^{kj_^_l ~ kj)z - lAt) + [v + lu] Apj_i exp [+i{kj_;^ - k'^)z + lAt) } 

- (fc, - fco) Ap,j (3.3.8) 



IV. COMPARISONS 



See "Wideband pulse propagation: single-field and multi- field approaches to Raman interactions" by Kinsler and 
New Q. 



V. SUMMARY 



I describe how to model a multi-frequency field such as that seen in a wideband Raman generation experiment 
using a single field envelope rather than a set of envelopes, one at each Stokes or anti-Stokes frequency. This requires 
that the field be propagated taking into account wideband effectSjas described by either the SEWA theory of Brabec 
and Krausz , or the more general GFEA of Kinsler and New . 

The usefulness of this single-field approach is not restricted to the Raman interaction described in this paper. It 
would be equally valuable for a near-degenerate optical parametric oscillator, or indeed any system where any two or 
more field components contain spectra that start to overlap as the pump or probe pulses get shorter. 

It is important to note that it will usually only be more efficient to use a single-field simulation if pump pulses are 
very short, and effects like the next-nearest neighbour field interactions, neglected in the multi-field theory, need to be 
included, or if the extra computational overhead is not inconvenient. This is because in a single-field simulation, a very 
fine time-resolution is necessary to model the polarization and field oscillations closely enough to get good numerical 
convergence. However, this situation improves when the Raman transition has a smaller frequency compared to the 
pump pulse frequencies. One useful side effect of the fine time resolution is that it naturally gives a wide spectral 
bandwidth, so that many Stokes and anti-Stokes lines are modeled quite naturally. Further, our single-field model 
could be invaluable in modeling a short pulse pump-probe experiment where the probe frequency does not match 
any of the Stokes or anti-Stokes spectral lines generated by the pump pulse(s). A multi-field simulation would then 
need arrays for both the pump and probe Raman 'ladders" of Stokes/anti-Stokes lines, and the role of next nearest 
neighbour interactions (ignored in the multi-field model) could well become more significant. 
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In summary, the advantages of our single-field approach are twofold. First, it includes more physics than the multi- 
field approach, even compared to a multi-field approach enhanced by adding GFEA corrections to the propagation of 
the field components. Secondly, it deals effortlessly with the complications of overlapping spectra in the multi-field 
case. 
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